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We use the Percus-Yevick approach in the chemical-potential route to evaluate the equation of 
state of hard hyperspheres in five dimensions. The evaluation requires the derivation of an analytical 
expression for the contact value of the pair distribution function between particles of the bulk fluid 
and a solute particle with arbitrary size. The equation of state is compared with those obtained 
from the conventional virial and compressibility thermodynamic routes and the associated virial 
coefficients are computed. The pressure calculated from all routes is exact up to third density 
order, but it deviates with respect to simulation data as density increases, the compressibility and 
the chemical-potential routes exhibiting smaller deviations than the virial route. Accurate linear 
interpolations between the compressibility route and either the chemical-potential route or the virial 
one are constructed. 
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I. INTRODUCTION 

The thermodynamic and equilibrium statistical- 
mechanical behavior of liquids and dense gases can be 
conveniently described by means of distribution func¬ 
tions [H-Q . The most widely used one is the pair correla¬ 
tion function or radial distribution function (RDF) g{r), 
which is particularly appropriate for the study of systems 
of particles interacting with a pairwise additive potential. 
For such fluids, there exists a set of well-established, rig¬ 
orous relationships connecting thermodynamic quantities 
with configuration integrals over RDF. In particular, the 
compressibility, energy, and pressure (or virial) equations 
provide well-known routes to the thermodynamic prop¬ 
erties of the fluid CHi. 

Thermodynamic properties can also be obtained from 
another route that connects g{r) with the chemical po¬ 
tential g through the charging process of a test parti¬ 
cle in the fluid [ 3 , 0 , 0 . This represents the chemical- 
potential route {g route), which has remained almost un¬ 
explored until recent years, when it was used to obtain a 
new Percus-Yevick (PY) equation of state (EOS) for the 
hard-sphere system 0. Subsequently, this has been for¬ 
mally generalized to arbitrary multicomponent systems 
and applied to additive hard-sphere (AHS) mixtures 0. 
In addition, the EOS of sticky hard spheres (or Baxter 
model 0 ) was derived and the critical point associated 
with a liquid-gas transition was captured by PY results 
in this way |ll| . 

All these routes (virial, energy, compressibility, and 
/i) are formally exact and thermodynamically equivalent, 
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provided that the exact RDF is employed 0. Actually, 
however, only approximate evaluations of g{r) are avail¬ 
able for systems in dimensions d > 1 with nontrivial in¬ 
teractions. In general, the RDF is obtained with some 
accuracy from numerical simulations 0, by solving in¬ 
tegral equation approximations (e.g., the PY [l^ and 
hypernetted chain approaches (la . fl4l| ) or from density 
functional theories mM- In this context, the g route 
has demonstrated to provide an alternative and useful 
path for the study of the thermodynamic properties of 
fluids. 

With the aim of extending the use of the g route and to 
gain insight into its properties, in this paper we apply the 
g method to find the EOS of the hard hypersphere fluid 
at spatial dimension d = 5. Besides the intrinsically in¬ 
teresting properties of hard particle systems, they are an 
important basis for constructing more complicate mod¬ 
els, so that there are active theoretical efforts to study 
them in dimensions d > 3 (see, for instance, Refs. [l7H^ 
and references therein). 

The g route is based on a charging process where a 
test particle (solute) with tunable interaction is inserted 
into the bulk fluid (solvent particles). As a consequence, 
the application of the method requires the solute-solvent 
RDE of the corresponding binary system in the infinite 
dilution limit. To our knowledge, the only systematic 
theoretical method for the evaluation of the RDE of 
AHS fluid mixtures at dimensions higher than d = 3 is 
the so-called rational-function approximation (REA) [s^ . 
Its simplest implementation provides the solution of the 
Ornstein-Zernike relation coupled with the PY closure 
for AHS mixtures in odd dimensionalities [s^. There¬ 
fore, we adopt the REA technique at the PY level in the 
present study. 

This work is organized as follows. In Sec.|TT]we present 



2 


^ = 0 C = 0.2 C = 0.4 



FIG. 1. (Color online) Cartoon of the charging process of the 
solute from a point particle = 0) to a particle equivalent 
to any solvent particle (^ = 1). 


the basic formulation of thermodynamics routes for hy¬ 
persphere systems. Section m gives, within the PY the¬ 
ory, the five-dimensional RDF for a coupled particle with 
arbitrary size, a key quantity for evaluating the EOS in 
the fi route. The technical details are elaborated in the 
Appendixes. The results are presented in Sec. |lVl Fi¬ 
nally, in Sec. 13 we offer some conclusions. 


II. FRAMEWORK 


The EOS for a single component system of hard d- 
sphere particles can be expressed as a relationship be¬ 
tween the compressibility factor Z = P/pk^T (where P 
is the pressure, p is the average particle density, fee is the 
Boltzmann constant, and T is the temperature) and the 
packing fraction 77 = Vdpo'‘^, with a the diameter of the 
particles and Vd = (7r/4)‘^'^^/r(l -|- d/ 2 ) the volume of a 
d-dimensional sphere of unit diameter. In the p route, 
the EOS is given by [sl, H, d, 

Z^( 77 ) =+ (2^^ _ 1)/ 


9 iv)-[ dttg{r]t) 
Jo 


( 2 . 1 ) 


with 


g(d)= 2dli _^ + (2.2) 


where g^rj- ^) = 312(0'12) is the contact value of the solute- 
solvent RDF, which depends on both 77 and the coupling 
parameter The latter is defined as the ratio between 
the solute and the solvent diameters, so that the mini¬ 
mum possible distance between solute and solvent parti¬ 
cles is 


1 + ^ 

ai2 = 


( 2 . 3 ) 


Thus, ^ regulates the strength of the interaction between 
the test particle and the rest of the fluid. When ^ = 


0 the test particle is a point that cannot penetrate the 
solvent particles, while when ■C = 1 the test particle is 
indistinguishable from any particle of the bulk fluid. This 
charging process is schematically illustrated by Fig. [T] 
For comparison, the compressibility factor in the virial 
and compressibility routes are expressed by 

Zv(77) = l + 2''-i773(77;l), (2.4) 

and 

Zc{v)= [ dtx~^{r]t), (2.5) 

Jo 

where xiv) is the isothermal susceptibility. Since hard 
spheres are athermal, the energy route to the EOS be¬ 
comes useless. 

In order to evaluate 312(0')*^) within the RFA approach 
[3^, it is convenient to introduce the Laplace functional 
defined by 


Gij{s)= / drr3y(r)6'„(sr)e (2.6) 

do 

where gij{r) is the RDF of the pair (7,3) and On(sr) is 
the reverse Bessel polynomial of order n = (d — 3)/2 [^ . 
This functional is directly related to the static structure 
factors (k) of a multicomponent fluid, 

S^j{k) = XiSij + 

where k is the wavenumber, Xi is the mole fraction of 
species f, and i is the imaginary unit. The functional 
Gij{s) provides us with all the necessary information 
about the structure and thermodynamics of the fluid 
state. In particular, the contact values are [s^ 

= lim (2.8) 

S—¥00 

Gij being the contact distance of the pair (7,3). 


III. PY APPROACH 


The exact solution of the PY approximation for d-odd 
AHS mixtures with any number of components has been 
described in detail in Ref. [s^ within the framework of 
the RFA method. For the sake of consistency, the main 
expressions particularized to Five-dimensional ( 5 D) bi¬ 
nary mixtures are presented in Appendix |3 From them 
one can take the limit of solute infinite dilution to obtain 
the solute-solvent contact value 3(77;^). The details are 
worked out in Appendix |B] and here we only quote the 
final result: 


3(3; 0 


2 + 377 - c + 3(1 - 77 + ce)c + - 3) 

(1+0^(2 + 377- 0(1 - 3 + CO 


( 3 . 1 ) 
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FIG. 2. (Color online) Plot of the solute-solvent contact value 
5 (^ 1 0 versus the solute-to-solvent size ratio ^ at a packing 
fraction 77 = 0.2. The values of and its slope at ^ = 0 

are exact, but the full curve is the approximate prediction of 
the PY theory [see Eq. mi]- 


where 


From here one can obtain the compressibility route to 
the EOS via Eq. (12.51) as 

^ 2(^-2 + 13577 -h I23O772 -h 3645773 -h 990 ? 7 '‘ -b 252773 

“ 22577(1 - 77)5 

(3.6) 


IV. RESULTS 
A. Virial coefficients 

We start this section by considering the virial expan¬ 
sion 

00 

= 1 + ^ 6„77”-3 (4.1) 

n=2 

for the 5D fluid as predicted by the three routes in the 
PY approximation. The virial coefficients bn'^ and bn^ 
corresponding to the virial and compressibility routes are 
obtained from Eqs. (13.311 [combined with Eq. (12.41) ] and 
(13.61) . respectively. As for the /r route, Eq. (12.11) implies 
that 


C = + (3.2) 

In the special case f = 1 one recovers the solvent-solvent 
value, namely. 




(^3 — 1 -b 3377 -b 8777 ^ -b 6773 
6677(1 — 77)3 


(3.3) 


In the opposite limit = 0 one obtains the exact results 
0 9{m^) = 1/(1 - 9 ) and %(7?;= 5??/(l - 
77 )^. Eigure [2] shows the ^-dependence of g{r];^) for the 
representative packing fraction 77 = 0 . 2 . 

Once ( 7 ( 77 ;^) has been analytically determined, the in¬ 
tegral in Eq. (12.21) can be evaluated with the explicit re¬ 
sult 


giv) = 


2 + St] 


C(5C-2) + 2 


62C3(2-b 377 -C) 11 -^ 

-5C + fc^ + ^C" + §C" + 2(i-C)2 


2 + 3t] 
1-77 


-c-rMin(i + c) 


(3.4) 


where ^ = ^/(l — 77 ). The compressibility factor in the /i 
route can be easily evaluated from Eq. dUD by numerical 
integration. 

Equation (13.311 can be used to evaluate the virial route 
to the EOS from Eq. (12.41) . Einally, the isothermal sus¬ 
ceptibility of 5D hyperspheres in the PY approach is 

[11 ill 


xiv) 


e 


[5(l-b 677 -b 3772)-2(2-b 377)C]. 


(3.5) 


where the coefficients which are defined by the ex¬ 
pansion 


gig) = 

n—0 


(4.3) 


are easily obtained from Eq. (EH). 

Table HI contains the first 12 virial coefficients obtained 
from the three routes in the PY theory. As is well 
known, the PY theory yields the exact virial coefficients 
up to third order, regardless of the thermodynamical 
route. However, discrepancies among results from dif¬ 
ferent routes appear for upper order coefficients {n > 4). 
As a peculiar feature, it may be observed that, in con¬ 
trast to the three-dimensional (3D) case 0 , the route 
yields irrational virial coefficients (due to the logarithmic 
term) for n > 5. On the other hand, the virial and com¬ 
pressibility routes yield rational numbers for the virial 
coefficients of hypersphere systems in all odd dimensions 


Table El shows the three sets of PY virial coefficients 
in numerical format and compares them with the exact 
result for n = 4 12311 and with recent accurate values 
for 5 < n < 12 [^. The relative deviations are dis¬ 
played in Fig. 121 It is quite clear that the compress¬ 
ibility route gives the closest agreement with the exact 
values, while the values calculated from the g and virial 
routes are rather similar, with a slight improvement of 
bif^ over bn\ It is interesting to remark that the three 
PY routes capture the alternating sign change between 
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TABLE I. First 12 virial coefficients obtained from the PY theory according to the virial compressibility (bn^), and ^ 

routes. 


n 


l(v) 




b 


M 


2 


16 


16 


16 


3 


106 


4 

5 

6 
7 


9 

10 
11 
12 


196 

1697 

2 

-2999 

164 989 
4 

-466 319 

184 953 797 
32 

1193 474 849 
16 

63 809 313 739 
64 

109 613124 547 
8 


106 

1459 

4 

2147 

2 

12 233 
8 

10 591 

3 790 541 
64 

24 293 447 
32 

1 124 833117 
128 

3 463 610 957 
32 

704 051293 633 
512 


106 

3561 

16 


2000 In 2 - 


1607 


303125, „ 19 465 513 

- 

47 221 875, „ 64427 077 

14 ~ 28 

2 998 209 375 , „ 49 556 094 191 

■ In 2 + - 


32 

2 345 096 875 In 2 - 


768 

87495 533 291 


54 

1756 815 571 125, 19 449143196 527 

-32-- 

431 924 172 778 125 224 301 227 394 011 

352 264 

6 786 922 855 865 625, „ 2 369 375 540 853 543 359 

■ In 2 H- 


256 


129 024 


TABLE 11. Comparison between the numerical values of the 
virial coefficients for 4 < n < 12 as obtained from several 
routes in the PY approximation and the exact values (b^^^) 
for n = 4 and 5 < n < 12 [33|. 


n 


10*bi''^/b”-^ 


lo^fo^rV^r' 

4 

759.7248 

478.5156 

890.5029 

543.3655 

5 

129.5219(16) 

129.4708 

163.8031 

129.7955 

6 

9.8184(19) 

-28.6007 

14.5829 

-23.3475 

7 

4.165(2) 

24.585 

6.313 

22.061 

8 

-1.127(4) 

-17.372 

-2.206 

-15.557 

9 

0.789(5) 

13.457 

1.768 

12.130 

10 

-0.468(10) 

-10.855 

-1.279 

-9.824 

11 

0.309(11) 

9.068 

0.984 

8.235 

12 

-0.23(2) 

-7.79 

-0.78 

-7.09 


n = 7 and n = 12, while a negative sign of be is wrongly 
anticipated by the virial and /i routes. In the PY case, 
the alternating character is related to the existence of 
a branch point singularity on the negative real axis (at 
r] = —3/2 + 5-\/3/6 ~ —0.0566), which determines the 
radius of convergence of the virial series [^ . 

Regarding the closeness between and bn '^, it turns 
out to be higher in five dimensions than in three dimen¬ 
sions. While in the 3D case the ratio bn'^/h^'^ mono- 
tonically decreases from 0.955 (n = 4) to | (n —oo) 
0, in the 5D case the ratio first increases from 0.881 
(n = 4) to a maximum value 1.225 (n = 6) and then 
(from n = 8) monotonically decreases towards an asymp- 



FIG. 3. (Color online) Normalized differences between the 
virial coefficients in the virial (squares), compressibility (cir- 
cles j m d p (triangles) routes, with respect to the exact values 

totic value 1.0646. It can then be speculated that the 
general similarity between and b^^ tends to increase 
with increasing dimensionality. We will return to this 
point at the end of Sec. IIVBI 
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FIG. 4. (Color online) Compressibility factors evaluated from 
computer simulations (symbols) and from an accurate 
Fade approximant ^[ 4 , 5 ] (dash-dotted line) [30l . 1^ are com¬ 
pared with PY results in the virial (long dashed line), com¬ 
pressibility (short dashed line), and /r (solid line) routes. The 
dotted vertical lines mark the freezing and melting packing 
fractions [^ . 


B. Equation of state at finite density 

As said before, the compressibility factor of 5D hy¬ 
perspheres in the PY /r route can be obtained through 
numerical integration from Eq. (EH) by making use of 
the explicit expression (13.41) . The virial and compress¬ 
ibility routes yield Eqs. ( 1 ^ . combined with Eq. da, 
and da, respectively. 

EigureS] shows the density dependence oi Z{t]) accord¬ 
ing to the PY theory in the virial compressibility 

(Z^^), and a (Z^'^) routes. Predictions of computer 
simulations are denoted by symbols, while the dash- 
dotted line represents a Z [4 5 ] Pade approximant pill 
based on hi gh p recision calculations of the first ten virial 
coefficients [^. The freezing (per® « 1.06, ry « 0.174) 
and melting (per® r; 1.25, rj Ri 0.206) densities calculated 
for the D 5 lattice are also indicated. As expected, all 
three PY routes converge to the exact results at very low 
density. On the other hand, as the packing is increased, 
the compressibility route overestimates Z, whereas the 
virial and /i routes underestimate it, the values from the 
/i route being slightly better than those from the virial 
route. 

The relative deviations of Z^'^, and Z^^ from 

the exact compressibility factor Zex are plotted in Pig. [S] 
It can be observed that the deviations of the compress¬ 
ibility and /i routes from computer simulations are very 
similar (but with different sign) within the fluid phase 
(p < 0.2). For instance, at p = 0.1, the deviations be¬ 



FIG. 5. (Golor online) Relative differences between the com¬ 
pressibility factor Zex evaluated from computer simulations 
(rj > 0.14) or from the Fade approximant Z[ 4 ^ 5 ] (p < 0.14) 
1^. 1^ and the FY results in several routes (filled symbols). 
Open symbols correspond to the interpolations Zv,c (squares) 
and Z^,c (triangles) from Eqs. (14.41) and (14.51) . with a = 0.44 
and Q = respectively. 

tween simulation or Pade-approximant results and the 
PY theory become 3.2%, 2.2%, and 2.5% for the virial, 
compressibility, and p routes, respectively, while they in¬ 
crease to 11.7%, 9.8%, and 9.3%, respectively, when the 
packing fraction p = 0.2 is reached. 

Looking for a better agreement with numerical ex¬ 
periments from PY solutions, one can construct new 
EOSs through Carnahan-Starling-like interpolations of 
the form 

Zv,e(p) = aZ^il) + (1 - a)Zr{v). (4.4) 

Z,,M = + (1 - a)Zr{v)- (4.5) 

To test this possibility, one can define density-dependent 
weight functions av.c(??) = - ZeAv)]/[Zf'^{r]) - 

Zv^iv)] and a^,c(p) = [Zf'^iv) - Zex(p)]/[2c^^(p) - 
Zj^^iv)] to check to which extent they are weakly depen¬ 
dent on rj. This is done in Fig. [51 which shows evaluations 
of the parameter a(p) from available data. Variations ob¬ 
served at the lowest densities (p < 0.07) are due to the 
ratio of very small values; actually, all three routes yield 
accurate results for very dilute gas conditions. On the 
other hand, it is clear from Fig. [5] than there is not a 
unique value of av,c(p) or in tho fluid phase re¬ 

gion (p < 0.2). However, it is to be noted that a = ^ 
and a = 0.44 may be reasonable choices for Z^_c and 
Zv,c, respectively. It can be observed from Fig. [5] that 
the value a = 0.44 in Eq. (14.41) is better in the region 
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FIG. 6 . (Color online) Ratio a = - Z^^jiZT^ - Z) 

with Z = (squares) and Z = Z^^ (triangles), using the 
compressibility factor Zex obtained either by simulations JM 
(filled symbols) or from the Fade approximant ^[ 4 , 5 ] [^. l40l| . 
As reference, horizontal lines denote the values a = a = 
0.44, and a = |, and vertical lines refer to the freezing and 
melting packing fractions 


M i 0.2 than the simpler value a = | proposed in Ref. 

. Data obtained from the hybrid compressibility fac¬ 
tors Zy^c (with a = 0.44) and (with a = i) are also 
included in Fig. [5] Both EOSs have discrepancies with 
respect to simulation results lower than 0.3% for rj < 0.2. 

The improvement of the /i route over the virial route 
observed in the virial coefficients (see Table El and Fig. 
ED as well as in the EOS at finite densities (see Figs, 
m and [5]) can be rationalized by the following heuristic 
argument Q . Comparison between the exact statistical- 
mechanical formulas (EH and (12.41) shows that, while the 
fj. route is related to the (weighted) average 3 ( 77 ) of the 
contact value g{r];^) in the range 0 < ^ < 1 , the virial 
route is directly related to the local value g{ri; 1). Since 
both g{ri\i) and its first derivative dg{r]-,^)/d^ are given 
exactly by the PY approximation at ^ = 0 (see Fig. 12 ), 
it seems reasonable that the average value 5 ( 77 ) is better 
estimated than the end point 3 ( 77 ; 1) by the PY approx¬ 
imation. In this respect, it is also worth noting that the 
weight function (1 -b in Eq. (12.21) increases with ^ 
(except in the one-dimensional case, d = 1 ), so that the 
average value g{r]) is more influenced by the approximate 
values of 5 ( 77 ; near f = 1 than by the quasiexact values 
near ^ = 0. This bias strongly increases with increasing 
dimensionality [in fact, (1 -b ^)‘^“^ —)► {2'^/d)S{^ — 1 ) in 
the limit d —> c»], which explains why the discrepancy 
between the PY 77 and virial routes is much smaller with 
d = 5 than with d = 3. 


V. FINAL REMARKS 


In this paper, we have considered an application of the 
Kirkwood coupling parameter method (or 77 route) to the 
thermodynamics of classical fluids. Specifically, the EOS 
of the 5D-sphere fluid has been derived from the 77 route 
by exploiting the knowledge of the exact solution of the 
PY theory for this system. The basic quantity required 
for this implementation is the contact value of the RDF 
for a partially coupled particle, which corresponds to a 
binary mixture with one infinitely dilute component. Al¬ 
though the PY approximation had been solved for the 
5D-sphere system by Freasier and Isbister [s^ and by 
Leutheusser [4^, their methods do not provide the RDF 
of a binary mixture required here. For this purpose we 
have used the RFA method which exactly solves the 
PY equation for hypersphere mixtures in odd dimension¬ 
alities. The RDF obtained in this way has an analytical 
expression and the resulting EOS in the 77 route is derived 
from a numerical integration over density. 

By examining the first few coefficients in the density 
expansion for the EOS and comparing them with avail¬ 
able exact or highly accurate values, it is seen that, al¬ 
though the PY theory is exact up to the third coefficient, 
it rapidly worsens for increasing density order. Specifi¬ 
cally, the compressibility equation gives in general the 
most satisfactory values, followed by the 77 route, which 
becomes somewhat better than the virial one. Contrary 
to the case of 3D-sphere systems Q , linear interpolations 
of the PY EOS from different routes do not improve the 
evaluation of the virial coefficients. 

Regarding the behavior at finite densities, the pres¬ 
sure predicted by the 77 route is observed to lie below 
to simulation data within the fluid phase. Interestingly, 
the deviations of the compressibility and 77 routes with 
respect to the numerical experiments are roughly of the 
same magnitude (but of opposite signs) in this region and 
both are lower than those from the virial route. An ac¬ 
curate fit to simulation data in the fluid regime (77 < 0 . 2 ) 
is obtained if the compressibility and 77 routes are added 
together with equal weights. An alternative accurate fit 
to numerical experiments may be derived by a similar in¬ 
terpolation between the virial and compressibility routes 
with appropriate weights. These results show a rather 
regular behavior of the three PY routes around the exact 
EOS within the whole fluid region. 
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Appendix A: RDF for a general five-dimensional 
binary mixture 

For an AHS binary mixture at c? = 5, the exact PY 
solution of the functional (12.61) may be expressed by [s^ 

= —^2—[(Lo + Lis + L2S^) • B , (Al) 

where Lm {rn = 0,1, 2) and B(s) are 2x2 matrices. More 
specifically, 

B(s) = \ + p [cDo(s). Lo + 0i(s) • Li + 02(s) • L2], (A2) 

where I is the unit matrix and 0m(s) (jn = 0,1,2) are 
diagonal matrices with elements 




where 


0m (^) — 






1=0 


J'- 


(A3) 


(A4) 


Furthermore, 


'1 1' 

77 

■IO77MI/MS , Mr 

eri (712 

1 1 

Ms 

(1 — 677)2 ' 1 — 677_ 

(712 (72 


crt crt 


2r]Me/M5 
1 — 677 

3 

(Jl 0’ia'2Cri2 
0'lO’2Cri2 


-2 -^2 1 3^/2 






1 — 677 

(1 + 24r7)M6 


4 ^ 4(1 - 677) 2 


CTi Cri(T2 
(Ti(T2 (t| 


cri2 
C^12 <^2 


(AlO) 


Qi = 


l0r]^Me/Mi 
1 — 677 
, 5?7/2A/5 


XlCTi X2Cr2 

Xicrf X2cr2 

XlCTi X20'2Cri 

^5^ ^ 


1-677 \xi(Jla2 X2 <jI \ 3M5 


CTl 0 
0 (72 

277 


Xl(Ji X2<J2 

Xiuf X2<J2 

' 4 4' 

Xl(Ji X20'2 
4 4 

Xl<Ji X20'2 


10rjMl/M5 Mr 


+ 


lOrj^ 

M| [ (1 — 677)2 
25T]^Me/Mi 
(1 — 677)2 


1 — 677 


Xiafa2 072(72 


1 — 677 


271(7® a;2(7|(7i 

Xiafa2 272(7® 


(All) 


Q2 = 7: 


1 

'1 o' 

377 

Xldf 272(7f 

3 

0 1 

Ms 

Xidf X2(J2 


20r]^MQ/M, 


1 — 677 


Ln — 3 


ri 1] 



r 4 

XidJ 

41 

X2(J2 

577/Ms 

27iCTf 

272(7|cti 

1 1 

, (A5) 

X 

4 

xiaf 

4 

X2(J2 

1-677 

Xiafa2 

X2(72 


Li =3 


e^i eri2 

1 

f2Me 

1 1' 

0’12 er2 

' 1-677 

K Ms 

1 1 


+ 3 


(7l (72 

(7i (72 


10 

A4 


4 4' 

XlCTf X2(J2 
4 4 

XiCTi X2<72 


L2 . (A6) 


Here, M„ = Xia^ + X2<J2 is the Trth diameter moment 
and 77 = v^pMr, is the total packing fraction. As for the 
matrix L2, it obeys the quadratic equation 


(A12) 

In general, the solution to Eq. (TATI) must be obtained 
numerically [35| . However, an analytical form is possible 
for a binary mixture with an infinitely diluted component 
(the solute), as shown in Appendix I bI 

The asymptotic behavior of the functional Gij{s) re¬ 
quired in the evaluation of the RDF contact value [Eq. 
(EH)] can be derived from m as 


Qo -f Qi ■ L2 + Q2 • L2 • (Po + Pi • L2) = 0, (A7) = hm = [L2 • B^^^] 


(A13) 


where 


Pn = 77 
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where we have called Bqo = lims_>.oo 
in Ref. IsHl that 


Boo — (C 5 • Lq + C 4 • Li 


where 


Cm = (-1) 


m+1 
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Appendix B: Test particle limit 

The expressions in Appendix [^ apply for any five¬ 
dimensional AHS binary mixture in the PY approxima¬ 
tion. Now, without loss of generality, we adopt the con¬ 
vention 7 = 1 and 7 = 2 for solvent and solute particles. 
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respectively, and denote by cti = cr and a 2 = the re¬ 
spective diameters. Next, we take the test particle limit 
for the solute, i.e., cci —>■ 1 , 0:2 —>■ 0 , so that —>• cr”. 

This drastically simplifies the matrices (IA8I) - (IA12I) . 

Let us write the matrix L 2 as 


L 2 — 


All Ai2 
A 21 A 22 


(Bl) 


C22 = 377(3 -I- 5^ -I- 277 )Aii -I- (1 - 6 r 7 )A 2 i. 

It can be easily checked that if the solute is equivalent 
to a solvent particle (i.e., ^ = 1 ), then the elements A 12 , 
A 21 , and A 22 coincide with An, as expected. 

Now we determine the contact values. Taking into ac¬ 
count that 77 = 7r^pcr®/60 and using Eqs. (IA5I) . (IA 6 I) . 
(IA15I) . and (IBID , one obtains 


The (1,1) element of Eq. (IA7I1 yields a closed quadratic 
equation for An. The physical solution is identified by 
the condition lim^_^o An = finite with the result 


1 - 1-477 (1-^)(1 - 677) 

2 O 77 2 (h^C 
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- 3 O 77 
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An Ai2 
0 0 

-y;(l-t-10Aii) 277 ( 1 - 1 -10 A 12 ) 

1-677 2(1-677) 

0 0 


, (Bll) 
■ (B12) 


Insertion of this into Eq. (IA14I) yields 


where C is defined by Eq. (IQ) . Since An is a coefficient 
related to the bulk fluid, it does not depend on ^ and 
its expression is equivalent to that of the one-component 
fluid [as given by Eq. (ElO) of Ref. [l^]. As for the 
remaining elements of the matrix L 2 , they are given from 
Eq. dSI by 


Q12/4C"- 277 ( 2 + 3e)An 
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Its inverse matrix is 
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Finally, from (IA13I) one has 
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5111^1 ) ,^2/y^^ _ 2077(2 + 377) 


(B15) 
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8 ( 1 - 677 ) 2 ^ 


where 


2077(2 -h 377 )Ai2 - 377(3 + 2r, + 50 
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12 
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ai2 =(1 + 0 " + 377(5 + 8^ - 2 O - 40) 

-h 1277^6-f 4^ - 60 + 30) + 1277^ (B6) 


32i(c7^2) = 


1 — 677 A2 


(1 + 0 /An - 2077(2 + 377) An 


, (B17) 


621 = 13-h30C-f 15O + 6?7(19-f20C-15O)+48?7^ (B7) 

022 =80 + 377(1 -f 8 ^ -h 180 - 20 O) 

-f 1277^(5-f 8^-120+ 60) +1277^ (B8) 

622 = 2-hl2C-hl5O + 677(6-hC-15O) + 67702-3O, (B9) 


322 (ctJ) = 


1 [ 2077(2 -I- 377 )Ai 2 — 3 ? 7(3 + 277 -|- 50 A21 


0 /Aii - 2077(2 + 377) 


A 22 


All 

(BI 8 ) 


As expected, gi 2 (< 7 /l) = 521 ( 0 '//). By inserting the ex¬ 
plicit expressions of An and A 12 it is straightforward to 
check that Eqs. (IB15I) and (IB16I) reduce to Eq. (13.31) and 
dS+D, respectively. 
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